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ABSTRACT: 


The optimization of nutritional parameters and concentrations for the protease production by Bacillus sp. in submerged fermentation 
was carried out using response surface methodology (RSM) based on the central composite design (CCD). The design contains a 
total of 20 experimental trials containing starch, soybean meal and CaCO; as model factors. The mutual interaction between these 
variables resulted into 1.48 fold increase in protease activity as compared to the mean observed response at zero level of all 


variables. 
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INTRODUCTION 

Proteases are the hydrolytic enzymes that cleaved 
protein molecules to release amino acids. Proteases 
have variety of sources like pants, animal, bacterial, 
fungal and viral origins. Protease production from 
Bacillus sp. was achieved by many researchers [1-3]. 
Protease has many applications in food, textile, paper 
and pulp, pharmaceutical, baking and beverages 
detergent and leather industries [4-5]. Industrially 
important enzymes including proteases have 
traditionally been obtained from submerged cultures 
because of easy handling, greater control of 
environmental and nutritional factors. The most 
frequently used operation in biotechnology is to 
improve the fermentation conditions for maximizing 
cell density, high of desired metabolic product or 
enzyme levels in microbial system. This approach is 
time consuming and also ignores the combined 
interactions between physical as well as nutritional 
factors [6-7]. In contrast, RSM includes factorial 
design and regression analysis that helps in evaluating 
the effective factors and in building models to study 
interaction and select optimum conditions of variables 
for a desirable response [8-9]. Recently, a number of 
statistical experimental designs with response surface 
methodology have been employed for optimizing 
enzyme production from microorganisms [10-11]. 


However, 3D and counter plots for response surfaces 
can provide a good way for visualizing the parameter 
interaction. Therefore, technique is often used for 
predicting optimum process conditions for microbial 
enzyme production. It is well known that extracellular 
enzyme production in microorganisms is greatly 
influenced by nutritional factors like carbon sources, 
nitrogen sources and mineral salts [12]. Enhancement 
in extracellular protease production from Bacillus sp. 
by optimization of fermentation conditions has not 
been attempted so far. Therefore, considering the 
many industrial applications of protease, we report 
here the optimization of extracellular protease 
production from Bacillus sp. as a result of the 


interactive effects of three variables (i.e. starch, 
soybean meal and CaCO3) using response surface 
methodology. 


MATERIALS AND METHODS 


Microorganism 

A strain Bacillus sp. a garden soil producing 
extracellular protease. The strain was maintained on 
2% nutrient agar slants at 4°C and as a glycerol stocks. 


Chemicals 

Chemicals and media used were of analytical grade 
purchased from the Fisher Scientific (Pittsburgh, PA 
USA). 


Composition of medium 

Starch- soybean meal- CaCO; (SSC) medium (pH 8.0) 
was used for the growth and production of protease by 
Bacillus sp. 


Inoculum preparation 

Seed inoculum was prepared by growing the isolate on 
Nutrient agar at 30°C for 24h. The cells were 
suspended in saline and cell density was measured 
spectrophotometrically (Shimatzu UV-2501 PC, Japan) 
at 600 nm. 


Experimental design and optimization by RSM 

The optimum levels for protease production by the 
Bacillus sp. in submerged fermentation with respect to 
starch, soybean meal and CaCO3, were obtained by 
single factor optimization. The experiments were 
conducted in 250 ml Erlenmeyer flasks containing 100 
ml of sterilized SSC medium inoculated with the 
freshly prepared 1% (10° cells/ml) bacterial suspension 
(as discussed earlier) and incubated for 12 hrs at 30°C 
under shaking culture condition (150 rpm). After 
fermentation, the cell-free supernatant was obtained by 
centrifugation at 10,000 rpm and the extracellular 
protease activity of the fermented broth determined. In 
the next stage RSM was used to study the interactive 
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effects of three variables, i.e. starch, soybean meal and 
CaCO, for improving total protease production. Each 
factor in the CCD was studied at three different levels 
(-1, 0, +1). The minimum and maximum ranges of 
variables were investigated and the full experimental 
plan with respect to their values in actual and coded 
form was listed in Table 1. 


Table 1. Experimental range and levels of the three 
independent variables used in RSM in terms of actual 
and coded factors 


Variables Range of levels 


A C A C A C 
X, - starch 10 -l 415 20 +41 


0 
%-SBM 10 -1 45 9 20 41 
0 


X3;—CaCO; 0.1 -1 0.3 05 +41 


A-actual values, C-coded values, actual values are expressed in g% 
SBM-Soybean meal 


A 2? factorial CCD proposed by Box er al.[13] with 
three factors leading to a total of 20 sets per 
experiment was formulated to optimize the starch, 
soybean meal and CaCO; concentrations. All the 
variables were taken at a central coded value 
considered as zero. The concentrations of these 
nutritional factors in the production medium were 
varied according to the experimental design as shown 
in Table 2. Experiments were conducted in duplicate. 
Using RSM, the relationship among the variables, i.e. 
starch, soybean meal and CaCO; were expressed 
mathematically in the form of a polynomial model, 
which gave the response as a function of relevant 
variables. The present work was based on the CCD to 
obtain the experimental data, which would fit an 
empirical, full second-order polynomial model 
representing the response surfaces over a relatively 
broad range of parameters. RSM had not only been 
used for optimization of medium components in the 
fermentation process [14] but also for studying the 
combined effects of culture parameters [15-16]. An 
empirical second-order polynomial model for three 
factors was in the following form: 
3 3 3 
Y>aQotsaixitS Seaijxixf on... {2} 


where, y was the predicted response (enzyme 
production) used as a dependent variable; xi (i = 1, 2 
and 3) were the input predictors or controlling 
variables; and do, ai (i = 1, 2, 3) and aij (@= 1, 2,3; j= 
i,..., 3) were the model coefficient parameters. The 
coefficient parameters were estimated by multiple 
linear regression analysis using the least-squares 
method. A second-order polynomial equation was then 
fitted to the data by least-squares optimization 
technique. This resulted in an empirical model that 
related the response measured to the independent 
variables of the experiment. 
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Assay of protease 
The cell free supernatant (1 ml) was mixed with 4 ml 
of casein (0.625% w/v) and incubated at 37°C for 30 
min. The reaction was stopped by addition of 5 ml of 
trichloroacetic acid (5%). Enzymatically hydrolyzed 
casein was measured by modified Folin Ciocalteu 
method [17], against casein treated with inactive 
enzyme as blank. A standard graph was generated 
using standard tyrosine solutions of 5 - 50 pg ml’. One 
unit of protease activity was defined as the amount of 
enzyme, which liberated 11g tyrosine per min at 37°C. 


RESULTS AND DISCUSSION 

A submerged culture was used for the production of 
extracellular protease enzyme from Bacillus sp. 
Preliminary experiments on protease production from 
the above strain indicated that the most important 
nutritional parameters were starch, soybean meal and 
CaCO3. Hence these three factors were considered as 
the independent variables and their effects on protease 
production were studied using a CCD of RSM. The 
results of CCD experiments for studying the effects of 
three independent variables, viz., starch, soybean meal 
and CaCOs, on protease production are presented in 
Table 3 along with the predicted and observed 
responses. The standard deviations on the observed 
responses are also presented in Table 3. 


Table 3. Observed responses and predicted values 


Protease activity (U ml! min’) 


Ru7=—SOti<“iS.... Residual Std. 
No. Observed response _ Predicted response Deviation 
1 86.00 82.9145 0.382 
2 93.00 103.6446 -1.318 
3 66.00 74.2737 -1.025 
4 122.00 117.8042 0.520 
05 178.00 183.8366 -0.0723 
6 197.00 190.5666 0.797 
7 187.00 191.6747 -0.579 
8 203.00 205.5056 -0.310 
9 109.00 97.3630 1.441 
10 149.00 142.4368 0.813 
11 141.00 140.6008 0.049 
12 144.00 141.8782 0.263 
13 187.00 175.9751 1.365 
14 195.00 199.4414 -0.550 
15 137.00 138.0140 -0.126 
16 136.00 138.0140 -0.249 
17 137.00 138.0140 -0.126 
18 136.00 138.0140 -0.249 
19 137.00 138.0140 -0.126 
20 136.00 138.0140 -0.249 
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Coefficients of the model given in eq. (1), were 
determined by the Gauss—Newton technique of least- 
squares optimization. 


The equation with the optimized coefficients is given 
by- 

y = 51.591 + 105.842 X; — 103.604X2 + 187.785X3 + 
22.800 Xj2 + 1,948X>3 = 35X73 = 34,.804X)) + 
23.989X 9 + 162.620X33 

where, y is the predicted response (protease activity), 
X, the starch in % X, the soybean meal in % and X; 
the CaCO; in %. The significance of the coefficient as 
determined by student’s t- test and p- values, which are 
listed in Table 4. The larger the magnitudes of t- value 
smaller the p- value, the more significant is the 
corresponding coefficient [18]. 


Table 4. Model coefficients estimated by multiple 
linear regressions 


Factor Coefficients Computed P- value 
t- value 
“Intercept = 51.591s««964— CSB 
X)- Starch 105.842 1.504 0.164 
X- SBM - 103.604 -1.457 0.176 
X3- CaCO; 187.785 1.961 0.078 
Xp 22.800 2.072 0.065 
X23 1.948 0.070 0.945 
X13 -35.000 -1.226 0.248 
Xu -34.804 -1.590 0.143 
X22 23.989 1.042 0.322 
X33 162.620 1.308 0.220 


SBM- soybean meal 


The results of the second order response surface model 
fitting in the form of analysis of variance (ANOVA) 
are given in Table 5. The fisher F-test with a very low 
probability value demonstrated a _ very high 
significance for the regression model [19]. The 
goodness of fit of the model was checked by the 
determination coefficient (R’). In this case, the value of 
the determination coefficient (R°= 0.988) indicates that 
only 1.2% of the total variations are not explained by 
the model. The values of the adjusted determination 
coefficient (Ad R= 0.956) is also very high, which 
indicates a high significance (p value < 0.01) of the 
model [20]. 


A higher value of the correlation coefficient (R=0.977) 
signifies an excellent correlation between the 
independent variables. At the same time a relatively 
lower value of the coefficient of variation (CV=4.4) 
indicates improved precision and reliability of the 
conducted experiments [18]. 
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The counter plots for response surface are plotted in 

figures 1-3 corresponding to the combined effects of 

soybean meal- starch, CaCO;- starch, CaCO3- soybean 

meal, respectively. The response surfaces obtained 

were suggesting that there was need of more starch 

concentration along with low soybean meal and CaCO; 
showing well-defined optimum nutritional conditions. 


Figures 1-3 revealed that all three substrates starch, 
soybean meal and CaCO; showed significant effect in 
protease activity. 


Table 5. Analysis of variance (ANOVA) for the three 
factorial design 
Variation SS DF MS FE p- 
Regression 27223.044 9 3024.783 46.381 0.000 
Residual 652.156 10 65.216 
Total 27875.200 19 


SS- sum of square, DF- degree of freedom, MS- Mean square, 
R’= 0.988, R= 0.977, Adjusted R= 0.956, * Significant at < 0.01 


Starch (%) 


10 06120~«WA 
Soybean meal (%) 


Fig. 1. Counter plot showing the mutual effect of soybean 
meal and starch on protease activity. CaCO; level was 0.3% 


Starch (%) 


025 0.30 
CaCO3 (%) 


Fig. 2. Counter plot showing the mutual effect of CaCO; and 
starch on protease activity. Soybean meal level was 2% 


Soybean meal (%) 


025° 030 035 0.40 
CaCO (%) 


Fig. 3. Counter plot showing the mutual effect of CaCO3and 
soybean meal on protease activity. Starch level was 2.0% 
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The optimum operating conditions obtained from the 
RSM models were starch of 2.0%, soybean meal of 2% 
and CaCO; of 0.5% with 203 Um! of observed 
protease activity and 205 Uml" of predicted activity. 
Total 1.48 fold increased in protease activity in both 
observed and predicted response in RSM as compared 
to central point at zero level, which corresponded to 
starch of 1.0%, soybean meal of 1.0% and CaCO; of 
0.3%. In recent years, there has been a great amount of 
research and development efforts focusing on the use 
of statistical methods, using different statistical 
software packages during process optimization studies, 
with the aim of obtaining high yields of protease in the 
fermentation medium [21-24]. The application of 
property designed approaches with multi-factor models 
allows process and biochemical engineers to design 
scale up strategies for increasing enzyme production. 


Thus, the RSM model performed well and offered 
stable response in predicting the combined interactions 
of the three independent variables, i.e. starch, soybean 
meal and CaCO; with respect to extracellular protease 
production. 


CONCLUSIONS 

This study compared the performance of the CCD 
using RSM in the estimation of fermentation 
performance for nutritional parameters and 
concentrations (starch, soybean meal, CaCO3) on 
extracellular protease production from Bacillus sp. 
Thus, RSM could be a very powerful and flexible tool 
for modeling the fermentation process due to 
corrective action arising from methodology and the 
associated estimation procedure. 
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